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Abstract 

We consider a new class of perturbation expansions, which incor- 
porate in a systematic way the available information about the diver- 
gent character of the perturbation series in QCD. The new expansion 
functions, which replace the powers of the coupling, are defined by the 
analytic continuation in the Borel plane, achieved through an optimal 
conformal mapping. We consider several possibilities of implementing 
the known behaviour of the Borel transform near the leading singulari- 
ties in the complex plane and using the corresponding expansions. These 
expansions have the same asymptotic convergence rate, but differ at low 
orders. We show that the new expansions allow a precise determination 
of the strong coupling as from the hadronic decays of the r lepton. 



1 INTRODUCTION 



As it is known, in 1952 Dyson obtained the famous result [T] that the renormal- 
ized perturbation series in Quantum Electrodynamics is divergent. During the 
subsequent decades, similar results have been obtained [2] -[9] for most of the 
physically interesting field theories, including Quantum Chromodynamcs, the 
modern theory of strong interactions. Dyson's proposal to regard a divergent 
series as asymptotic to the function sought for is nowadays generally adopted. 
This set a challenge for a radical reformulation of perturbation theory. The 
Feynman diagrams and renormalization theory yield, at least in principle, the 
values of all the expansion coefficients. They can tell us whether the series is 
convergent or not, but what we want to know is under what conditions the 
physical function can be determined from them. If the series were convergent, 
the knowledge of all the expansion coefficients would uniquely determine the 
expanded function. On the other hand, there are infinitely many functions 
having the same asymptotic expansion. A crucial task is to find effective addi- 
tional inputs that would be able to reduce or, if possible, remove the ambiguity. 
For instance, general properties such as causality and unitarity may be useful. 
For a recent discussion of the the ambiguity of field correlators represented by 
asymptotic perturbation expansions see [TO] . 

Along these lines, a possibility of exploiting better the low order expan- 
sion coefficients known at present from Feynman diagram is to combine them 
with the available information on the large-order behaviour of the expansion. 
Such an attempt was made in [llj . where a new expansion was defined by 
the analytic continuation in the Borel complex plane achieved by a conformal 
mapping^. Specifically, we use an optimal conformal mapping, as defined in 
[12j . As shown in [T3l [T4] , the new perturbation expansion separates the 
divergent character of the series from a part on the intrinsic ambiguity of the 
perturbation theory, which is generated by the infrared regions of the Feyn- 
man diagrams. This ambiguity is solved by choosing a prescription included 
in the definition of the expansion functions. 

In the present paper we discuss these ideas using for illustration the 
so-called Adler function [15]. This function is of interest in QCD for the de- 
termination of the strong coupling Og at a relatively low scale, the mass of the 
r lepton [l6]-[22]. In section [2] we briefly review the standard perturbation 
expansion of the Adler function in massles QCD. In section [3] we present the 
new perturbation expansion proposed in [11] and in section [5] we investigate 
a more general class of expansions functions, which have the same asymptotic 
convergence rate but differ at low orders. The convergence of the new expan- 
sions is checked using a realistic model proposed recently in the literature [20] . 

^For QCD, the idea of a conformal mapping in the Borel plane was suggested in [7] and 
was applied in several works (for references see [lljl. 
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In section [5] we briefly discuss the determination of as(m^), confirming the 
precise value reported recently in |21j . 



2 STANDARD PERTURBATION EXPANSION 

The Adler function [15] in massless QCD is the derivative 



where n(s) is the invariant amplitude of the correlator 

i [ d'^xe'i-^ < 0|T{y^(x)K(0)t}|0 >= {qf,q, - gf^.q^Ms). (2) 



Here s = is the energy squared and is the vector current for light quarks. 

In perturbative QCD, the reduced Adler function D{s) = D{s) — 1 is 
written formally as the renormalization-group improved series 

D{s) = J2Kniasis)r, (3) 

n>l 

where as{s) = as{s)/'K is the running coupling. Explicit calculations of Feyn- 
man diagrams in MS renormalization scheme with three flavours (nj = 3), 
predict the first four coefficients [23]-[25] 

Ki = 1, K2 = 1.639, = 6.371, = 49.076, (4) 

while an estimate of the next coefficient 

K, 283 (5) 

is available |20j . 

As discussed above, the series ([3]) is divergent, the coefficients displaying 
a factorial asymptotic increase, ~ n! |8j, |9j. In writing ([3|) we follow the 
convention often adopted in physical papers, writing the sign of equality even 
if the series on the right hand side is divergent and the equality is impossible. 
Analogous series in the present paper are understood in the same sense. 

The running coupling as{s) is obtained by setting the scale //^ = s in the 
renormalized coupling as{fM^) = as(/i^)/-7r, which satisfies the renormalization 
group equation 

^^d^ = /3(a.(^2)) = - 5:/3,(a.(M^)r+2 (6) 

The first two coefficients of the (3 function are universal (scheme-independent), 
and the next two were calculated in the MS scheme with n/ = 3 to four loops 
[261 127] . Therefore, the known coefficients are: 

/3o = 9/4, /3i = 4, /32 = 10.06, /Sg = 47.23. (7) 
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3 NEW PERTURBATION EXPANSIONS 



To define a new series, we start from the Borel transform B{u) of tlie Adler 
function, defined by the power series 

oo 

B{u) = h^u'^ , (8) 

n=0 

with hn defined in terms of the perturbative coefficients by: 

hn = ^, n>0. (9) 
Po n! 

According to present knowledge, the function B{u) has branch point singular- 
ities in the n-plane, along the negative axis - the ultraviolet (UV) renormalons 
- and along the positive axis - the infrared (IR) renormalons [El [9]. Specifically, 
the branch cuts are situated along the rays u < —1 and u>2. The nature of 
the first branch points was established in [6] and in [28] (see also [2Q]). Thus, 
near the first branch points, i.e. for u ~ —1 and m ~ 2, respectively, B{u) 
behaves as 

B(u) ~ 7 — B(u) ~ , , (10) 

^ ' (l+ti)Ti ^ ^ (l-n/2)T2' ^ ^ 

where the residues ri and r2 are not known, but the exponents 71 and 72 can 
be calculated using renormalization group invariance [HI [28l [20] , and have the 
values 

71 = 1.21, 72 = 2.58. (11) 
The series ([3]) can be formally written as the Borel-Laplace transform 

00 

D(s) = — I e-^/(l^oasis)) .-^2) 

/3o J 


Actually, due to the singularities of B{u) along the positive axis, the integral 
([12]) does not exist. The ambiguity in the choice of prescription is often used 
as a measure of the uncertainty of the calculations in perturbative QCD. It is 
convenient to define the integral by the Principal Value (PV) prescription: 

00 

b{s) = ^FY f e-"/Wi'^={^)) B{u) du , (13) 

Po J 


where 

00 00 
PV / f(u) du = lim - / [f(u + ie) du + f(u - ie) du] . 

J e-i-O 2 J 
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As discussed in |29] the PV prescription is the best choice if one wants to 
preserve as much as possible the analyticity properties of the correlators in 
the s-plane, which are connected with causality and unitarity. 

In order to define a new perturbative expansion of the Adler function 
we shall apply the method of conformal mappings |12) . This method is not 
applicable to the series because D{s) (regarded as a function of as{s)) is 
singular at the point of expansion as{s) = 0. The method can, on the other 
hand, be applied to ([8]), because B{u) is holomorphic at u = 0. 

We note that the expansion ([8]) converges only in the disk |u| < 1. A 
series with a larger domain of convergence can be obtained by expanding B{u) 
in powers of a new variable. As demonstrated in [12], the optimal variable, 
which leads to the best asymptotic convergence rate, coincides with the func- 
tion that performs the conformal mapping of the whole analyticity domain of 
the expanded function onto a disk in the new complex plane. 

It is generally assumed that B(u) is analytic in the w-plane cut along 
the real axis for u > 2 and u < —1 (for more remarks see [10^ I21j). Then the 
optimal variable defined in [12] reads [11] 



7r+^-yr^V2 

\/l + u + ^yl - u/2 

This function maps the n-plane cut for n > 2 and u < —1 onto the unit disk 
\w\ < 1 in the complex plane w = w{u), such that w{0) = 0, w(2) = 1 and 
w{—l) = —1. Then, the expansion 



B{u) = Y.dnw'' (15) 

n>0 



converges in the whole disk < 1 [12\ 111] . Moreover, as shown in [T2], the 
expansion (jlSp has the best asymptotic rate of convergence compared to all 
the expansions of the function B(u) in powers of other variables. 

The series (I15p can be used to define an alternative expansion of D{s). 
This is obtained formally by inserting (jlSp into (jl3p and interchanging the 
order of summation and integration. Thus, we adopt the modified expansion, 
defined as [II]- [H] 

D{s) = J2 dnWnis), (16) 
n>0 

in terms of the expansion functions 

oo 

Wn{s) = ^PV / e-"/(/^o'^^(^)) w'' du . (17) 
Po J 







We emphasize that the expansion (|15p exploits only the location of the 
leading singularities in the Borel plane. The series is expected to describe also 
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the nature of the singularities if a large number of terms is used. However, since 
the behaviour near the first singularities is known, cf. eq. (jlOp. it is convenient 
to incorporate it explicitly. This is achieved, for instance, by expanding the 
product (1 + w)^'^! (1 — B{u) in powers of the variable w: 

il + wf'''{l-wf"'^B{u) = Y,Cnw''- (18) 

n>0 

Here we took into account the fact that from the expression (jl4p of w = w{u) 
it follows that (1 — u/2) ~ (1 + w)'^ and (1 + ti) ~ (1 — w)'^ near the points 
w = —1 and w = 1, respectively. 

The expansion (jlSp converges in the whole disk \w\ < 1, i.e. in the 
whole cut complex n-plane. Moreover, since the singular behaviour of B{u) at 
the first branch points is compensated by the first factors in (jlSp . the series 
is expected to converge faster than (llSp . Also, the behaviour near the first 
singularities holds even for truncated expansions, which are used in practice. 
This suggests the definition of the new expansion 

D{s) = J2cnWn{s), (19) 

n>0 

where the expansion functions are 

oo 

with w = w{u) defined in ([HI. Actually, as discussed in [UlllT], the inclusion 
of the explicit behaviour near the leading singularities is not unique. The 
problem will be discussed briefiy in the next section. 

The expansions ([16]) and ([T9j) reproduce the coefficients Kn of the usual 
expansion ([3]) , when the functions ([T7j) and ([20]) are expanded in powers of the 
coupling. As shown in [T3], the new expansion functions are formally repre- 
sented by divergent series in powers of the coupling, much like the expanded 
correlator itself. 

A detailed comparison of the standard and the new expansions (jl6p and 
()19p was performed in [21], using as reference a realistic model for the Adler 
function proposed in [20]. This model, which expresses the Borel transform in 
terms of a few singularities, allows the calculation of the exact Adler function 
and of its perturbative approximants. To illustrate the comparison, we shall 
use the contour integral 

<5(0) = J_ I ^^^s)D{s), 



-v^.-v-., (21) 

\s\=so 
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Perturbative oder N 

Figure 1: Values of S^^^ for the model discussed in [20], obtained with the 
standard and the new expansions, as a function of the perturbative order A^. 
The horizontal band is the exact value from [20j with an experimental error. 
The calculations are performed with as(m^) = 0.34/7r. 



where sq = and uj[s) = 1 — 2s /sq + 2{s/sqY ~ {s/sq)^ is a kinematical 
factor. The quantity 5^^^ can be determined experimentally in terms of the 
total hadronic width of the r lepton [16]. Therefore, the integral (j2ip is of 
interest for the experimental determination of the strong coupling from the 
hadronic r-decays. 

To perform the comparison, we calculate 5'^^^ using the exact expression 
of the Adler function in the model [20] , and its expansions truncated at a finite 
order A^. The running coupling as{s) is obtained by integrating numerically 
the renormalization group equation ([6]) step by step along the circle |s| = m^, 
in terms of a reference value, taken to be as(m^). This procedure, which 
avoids the large imaginary logarithms appearing in expansions at a fixed scale 
is known as "contour improved" expansion of 5*^'^^ [16j. For illustration we use 
as{ml) = 0.34/7r, as in pi.^H]. 

The results are shown in Fig. [H where the gray band represents the 
exact result, to which we attached for convenience an uncertainty equal to the 
experimental error (see section [5]). As seen in Fig. [U the standard expansion 
is not able to reproduce the exact value at low truncation orders A^. Moreover, 
for higher the results start to oscillate wildly, due to the divergent character 
of the series. On the other hand, our new expansions (fT6]l and (fT9]) reproduce 
more closely the exact value. The expansion ()19p . which includes explicitly 
the singular behaviour at the first singularities, gives very good values even 
for low A'^, while the expansion (I16p gives a good approximation only at larger 
truncation orders A^. This is due to the fact that, in the model proposed 
in [20], the strength of the first singularities, expressed by the residues in 
pop , is quite important. Therefore, a large number of terms are required in 



' 1 ' 1 ' 1 ' 1 ' 1 ' 

□ standard expansion, Eq. (3) 
• new expansion, Eq. (16) 
■ new expansion, Eq. (19 ) 
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order to describe them, if the singularity is not exphcitly factorized. In other 
models, with milder leading singularities, the expansion (jl6p approximates 
well the exact results even at low values of N. Both expansions, (jl9p and 
(|16p . are actually much better than the standard expansion ([3]) in powers of 
the coupling. 



As remarked above (see also |1H l21j). while the expansion (|15p is unique, the 
explicit inclusion of the first singularities of B(u) contains some arbitrariness. 
The description of the singularities by multiplicative factors is a possibility, 
but is not a priori necessary. Moreover, the factors are not unique. For a large 
number of terms in the expansion the form of these factors is irrelevant, but at 
low orders one prescription may be better than another. In ()19p the dominant 
behaviour was included by simple singular factors expressed in the w variable. 
In this section we briefly discuss several other possible expansions. A more 
complete analysis will be presented elsewhere. 

To construct alternative approximants, we first remark that the dominant 
behaviour may be included by singular factors expressed in the u variable, 
instead of w. One may further multiply the leading factors by other functions 
analytic in the ii-complex plane cut along the real axis for u >2 and u < — 1. 
In particular, one may consider singularities on an unphysical Riemann sheet, 
or placed at u = 3 and u = —2, which are expected to occur as next-to-leading 
renormalons in the physical case. The additional factors can be expressed 
either in the variable u or in the variable w. We consider for illustration a set 
of new expansions written in the compact form 



4 OTHER ADMISSIBLE EXPANSIONS 





n>0 



where the expansion functions are defined as 



oo 



Wi(,) = -lpv/e 



w 



u/{Poas(s)) 



dn, 



(23) 







with the specific choices 




(24) 
(25) 
(26) 



(1 + wf^' (1 - wf^'' (2 - wf{3 - wf, 
(1 + wf^' (1 - wf^^ (1 + n/2)2(l - u/3f 



jJl V WkJ 



(27) 
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Figure 2: Values oi6^^^ for the model [20], obtained with various expansions, as 
a function of the perturbative order A^. The horizontal band is the exact value 
with an experimental error. The right panel enlarges the figure for N < 7. 



Mu) = {l + wf^'{l-wf^'l[ (l- — ) ' . (28) 

jj-^ V WkJ 

Here w = w{u), where the function w{u) is defined in and we denoted 
by wi = w{3 + ie), W2 = w{3 — ie), W3 = w{—2 + ie), W4 = ■w{—2 — ie) the 
positions the points u = 3 it ie and u = —2 it ie on the circle \w\ = 1 in the w- 
plane. In the absence of further information, the exponents of the additional 
singularities in F^{u), Ti{u) and F^{u) are arbitrary. 

In Fig. [2] we test the eficiency of the above approximants, using the 
model [20] and as{ml) = 0.34/7r as in the previous section. The results shown 
in the left panel confirm that at large all the choices become equivalent, 
as expecteqj- At low the differences between various approximations are 
more pronounced, but starting from N = 5, which corresponds to the present 
status in QCD, the range becomes rather narrow. In particular, for N = 5 
the expansion (I19p adopted in [21j gives the central value, the spread around 
it being ±0.006, close to the experimental uncertainty (we recall that for com- 
parison we attached an "experimental" error to the exact result, as shown by 
the gray band in Fig. [2]) . We shall discuss in the next section the implications 
of these results on the determination of a^- 

5 DETERMINATION OF a. 

The precise determination of Og from the hadronic decays of the r lepton is 
one of the most important results in perturbative QCD (for a recent review 

possible divergency at still higher values of A'^, not excluded in principle, is discussed 
in [T3[ll[2T]. 
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see [22] ) . The problem was revisited recently \19\ [20l [2T] , after the calculation 
of the Adler function to four loops [25| , the same order at which the f3 function 
of the renormalization group equation is known \26\ [27] . 

As mentioned in the previous section, the quantity relevant for the ex- 
traction of as(m^) is the integral (PT|) . We shall adopt the phenomenological 
value quoted in [20] : 



the renormalization group equation ^ along the integration contour. 

Taking the average over the various functional forms discussed in the 
previous sections, we obtain 



where the error includes the experimental uncertainty quoted in (j29p . the effect 
of a 50% variation of the coefficient K^, the uncertainty of renormalization 
scale [21], and an additional error of ±0.003 due to the freedom in choosing 
the expansion. 



The new perturbation expansion for QCD observables, proposed in [TT], is 
based on an optimal expansion of the Borel transform, which reproduces or- 
der by order the known perturbative coefficients calculated from Feynman 
diagrams and converges in the whole complex u-plane cut along two lines of 
the real axis. While the optimal expansion variable w is unique, a freedom 
exists in the way of including additional information, like the known behaviour 
of the Borel transform near the first branch points. The inclusion of this be- 
haviour is neither unique nor precisely formulated and we have no optimal 
method, in contradistinction to the choice of the optimal w, which is based on 
an "optimality" theorem [12j . 

In the present paper we investigated a set of possible ways of implement- 
ing the behaviour of the Borel transform at the leading singularities and tested 
their properties in the frame of a model proposed recently in [20]. The results 
show that, for the number of perturbative terms calculated at present from 
Feynman diagrams in QCD, the various approximants give very consistent 
results. 



Cen = 0-2042 ± 0.0050. (29) 





(30) 



6 CONCLUSIONS 
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Our analysis shows that the freedom in the implementation of the known 
singular behaviour of the Borel transform does not affect practically the extrac- 
tion of as from r hadronic decays: the central value of as(m^) is unmodified, 
the only effect being an additional error on of about ±0.003, comparable with 
the experimental error. Our final result is given in (j30|) . A detailed analysis, 
including a more general class of approximants and the investigation of alter- 
native physical models for the QCD correlator, will be presented in a future 
work. 
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